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ABSTRACT 

Type la Supernovae (SNe la) spectra are compared regarding the coefficient of the 
largest wavelet scale in their decomposition. Two distinct subgroups were identified 
and their occurrence is discussed in light of use of SNc la as cosmological probes. 
Apart from the group of normal SNe, another trend characterised by intrinsically 
redder colours is consisted of many different SN events that exhibit diverse proper- 
ties, including the interaction with the circumstcllar material, the existence of specific 
shell-structure in or surrounding the SN ejecta or super-Chandrasckhar mass progen- 
itors. Compared with the normal objects, these SNe may violate the standard width- 
luminosity correction, which could influence the cosmological results if they were all 
calibrated equally, since their fraction among SNe la is not negligible when perform- 
ing precision cosmology. Using largest wavelet scale coefficient in combination with 
long-baseline B — I colours, we show how to disentangle SN intrinsic colour from the 
part that corresponds to the reddening due to dust extinction in the host galaxy in 
the SALT2 colour parameter c, discussing how the intrinsic colour differences may ex- 
plain the different reddening laws for two subsamplcs. There are wavelength intervals 
for which the measured largest scale coefficient is invariant to the additional extinc- 
tion applied to a spectrum. Combination of wavelet coefficients measured in different 
wavelength intervals can be used to develop a technique that allows for estimation of 
extinction. 
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1 INTRODUCTION 

Type la Supernovae, widely accepted to be thermonuclear 
explosions of carbon/oxygen white dwarfs, represent a ho- 
mogeneous class of stellar objects, both spectroscopically 
and photometrically. The possibility to standardise the peak 
luminosity of SNe la, no matter how far they are, makes 
them efficient cosmological tool. Unfortunately, small differ- 
ences among SNe la directly translate onto errors on the 
cosmological parameters. We have reached a point where 
these errors are no more small for our goals. In order to 
obtain a better accuracy of distance determination, we aim 
to model the particular differences in a way to reduce as 
much as possible both systematic and statistical errors on 
the cosmological parameters. Indeed, several different tech- 
niques for fitting SN light-curves have been w i dely used (e.g. 
Jha. Riess fc Kirshnerll2007l ; iGuv et al. 1120071 ; IConlev et all 
20081'). Besides, some spectral indicators, such as line ra- 
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tips dNugent et al. Il995l) or pseudo-equivalent widths (EW ) 
l|Hachinger. Mazzali fc BenettillioOrj : bronder et al. ll200sT ). 



are found to correlate with absolute magnitudes, thus pro- 
viding a basis for an independent calibration method of 
the SN luminosity. A detailed study on the selection of 
global spectral indicators can be found in iBailev et ai~ll2009l . 
However, the spe ctroscopic diversity o f SNe la is multi- 



dimensional (e.g. iBenetti et al. 1 120051 : IWang et af~l 120091 : 
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iBranch. Dang fc Baron|[2009l ) . and a whole range of reported 
subclasses suggests that our understanding of models that 
correspond to these events is still not complete. Particu- 
larly, the observed spectral differences between SNe la can 
be linked to variations in their effective temperatures, which 
are associated with different amount of 56 Ni produced in the 
explo sions, i.e. different luminosities (see e.g. iMazzali et al. I 
l2007f ). 

In many cases the high-resolution SN spectroscopy is 
rare and, therefore, we cannot track the possible signature 
of the interaction between SN ejecta and surrounding cir- 
cumstellar material (CSM), normally exhibited by emis- 
sion lines in the spectra. It is demonstrated in this paper 
that the coefficient that corresponds to the largest scale 
in the wavelet decomposit ion of SN spectrum (for details 
see lArseniievic et al. I MPOST ), coupled with the SALT2 colour 
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(|Guv et al. 112007ft , is a parameter that provides additional 
information on SN subclassification, distinguishes SNe that 
show spectral peculiarities and most likely indicates different 
SN progenitor scenario and/or explosion mechanism. The 
wavelet inverse of the above-mentioned coefficient depicts 
the overall shape of the SN spectrum. Depending on the 
wavelength interval, the measured coefficient can be more 
sensitive to reddening. A method based on comparison of 
two largest scale wavelet coefficients from different wave- 
length intervals that indicates peculiar SNe, also those with 
anomalous extinction, is described. 

Both intrinsic SN colour and host galaxy dust extinc- 
tion are entangle d into single SALT2 colour parameter c 
(|Guv et al. 112007m . The unfolding of these effects which pro- 
vides a possibility to estimate the SN host galaxy reddening, 
E(B — l / )h ost , will be also discussed. 



2 DATA SET 

The discrete wavelet transform is applied to a sample of 73 
nearby SNe, mostly found in the literature and presented in 
Table [T] In addition, an unusually bright SN 2003fg (also 
known as SNLS-03D3bb. iHowell et al. 112006ft . observed at 
intermediate redshift z = 0.244, is also considered. All the 
low-resolution spectra (~600) have been deredshifted, but 
no reddening correction has been made. For most of the anal- 
ysis our attention is restricted to the rest-frame wavelength 
interval 3400 < A < 7000 A and epoch within [-10, +10] 
days relative to B-maximum. We perform then the Discrete 
Time Wavelet Transform (DTWT) of the spectra following 
the prescription given in lArseniievic et alTI 120081 . The de- 
composition is produced using Daubechies' extremal phase 
wavelets with 4 vanishing moments (D8). All SNe in our 
sample have relatively good photometry, and in order to ap- 
ply a consistent procedure, the SALT2 light-curve fitter is 
used for all supernovae to obtain the photometric parame- 
ters. However, the spectra employed in this study come from 
various sources and are not quite homogeneous, thus there 
are some limitations and uncertainties on spectra imposed 
by reduction issues that may affect the SN colour terms and 
certainly are not subject to intrinsic SN properties. 



3 METHOD AND RESULTS 

SN la spectra are decomposed into different wavelet scales 
following the procedure described in lArseniievic et"aL 1120081 . 
In order to normalise better the SN spectra, the authors do 
not include the contribution from scale 15, i.e. the coarsest 
wavelet scale in the spectral decomposition. Indeed, after 
measuring the featureless inverse of the wavelet coefficient 
from scale 15 at central wavelengths of B and V bands on 
our spectra, we see that a contribution from the largest scale 
to the B — V colour can be expressed as —0.0021 x wltcoeff 15, 
which gives a maximum amount of ~ 0.01 mag when applied 
to a range spanned by the coefficients, between —5 and (see 
Fig. [T(a7| ). The inverse of the wavelet from the largest scale 
may be regarded as a black-body curve or a sort of con- 
tinuum. On the other hand, we deal here with the spectral 
decomposition that does not map exactly to the underlying 
physics. In other words, there is no true physical meaning of 



each component (scale) in a wavelet decomposition; however 
one is able to measure the weight of each scale in terms of 
energy, as the wavelet power spectrum suggests. Further, a 
spectral analysis that is based on the wavelet decomposition 
certainly assures consistency of the method. Nevertheless, 
plotting the wavelet coefficient from the largest scale as a 
function of the SALT2 colour parameter, i.e. long baseline 
B — I colour in Fig. 1(b) leads to an interesting correlation. 



The robust regression fit in Fig. 1(a) reads c(SALT2) = 
(0.406 ± 0.025) + (0.125 ± 0.009) x wltcoeff 15, with the stan- 
dard deviation of fit a = 0.25, computed as the median 
absolute deviation of the residuals divided by the constant 
0.6745. In order to reduce large data scatter, two perceivably 
distinctive SN trends are identified, represented by solid and 



dashed lines in Fig. 1(a) 



c ( salt2) = (0.300 + 0.025) + (0.112 + 0.012) xwltcoeffi 5 (1) 
and 

c ( salt2) = (0.632 ±0.026) + (0.149 + 0.008) x wltcoeff 15 (2) 
respectively, with the robust measures of spread equal to 



0.09, i.e. 0.11. Similarly, in Fig. 1(b) the scale 15 coefficient 
is well correlated with a long-baseline B — I colour at B- 
band maximum corrected for Galactic reddening assuming 
E(B-I) = 2.35E(B-V). To calculate the uncertainty of the 
(B — I) values, the photometric errors with the uncertainties 
due to the fitting procedure are combined. 



A trend described by solid line in Fig. 1(a) is referred to 
as "normal" SNe. On the other hand, the earlier mentioned 
SN 2003fg is identified as a member of SN subclass (repre- 
sented by dashed line), that includes objects whose colour 
appears to be redder than the typical SN la, according to 
the fits from Eqns. [T] and [21 Looking closely to the members 
of this SN group, represented by filled circles, there is some 
doubt whether these differences reflect the different explo- 
sion scenario or distinct progenitor channels, whether it can 
be attributed to the existence of (overdense) shell-structure 
in or around the SN ejecta, the interaction with the circum- 
stellar material, massive progenitors, or finally, to different 
viewing angles of asymmetric envelopes, as will be discussed 
in subseque nt section. 

Indeed, iBlondin et al. 1 120061 found that spectrum of 
SN 2003fg matches best with the spectrum of SN 1989B 
at t = +3.5 days (with constraints on epoch or red- 
shift), an object also ide ntifi ed to belong to the sa me 
group. IHowell et al~ll2006l and iHillebrandt" et al. I l2007l de- 
bate between the possible explosion scenarios for SN 2003fg: 
a super-Chandrasekhar-mass progenitor explosion or the 
asymmetric explosion of a Chandrasekhar-mass (Mch) 
white dwarf? Som ething similar can be said for overlu- 
minous SN 2007if |Scalzo et al. I [2OI0I ) and for SN 2006gz 
(|Hicken et al. Il2~007ft , both with progenitor mass larger than 
Mch and both surrounded by an envel ope of unburned 
carbo n/oxygen. Interestingly, SN 2009dc (|Yamanaka et aL 
2009), another SN with a massive progenitor, seemqj 



to be closer to norm al SNe la in Fig. 1(a) Recently 
ISilverman et al.ll201ll published 9 spectra of this SN that 



1 A value close to 1.5 for the wavelet coefficient fr om the largest 
scale is measured from the available spectra from rTanaka et al. I 
120091 . 
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Figure 1. (a) Mean value of the wavelet coefficient from the largest scale for different spectra of nearby SNe la. The solid line corresponds 
to the best fit to "normal" SN data, while the dashed line describes the other SN class (see text for details), (b) Same as in (a), but using 
a long-baseline B — I colour instead of the SALT2 c parameter. Seven normal SNe that suffer negligible host extinction are identified as 
crossed symbols. 




cover the epochs from -7 days to +281 days and extend to 
3500 A. An additional check for classification of this ob- 
ject on the spectrum at -7 days shows that the scale 15 
coefficient for SN 2009dc is roughly —4 ± 1, therefore this 
SN belongs to the second group, as was expected. SNe la 
found to deviate from the main SN trend in Fig. [1] whose 
colours are described by Eqn. [2l may show significant dis- 
crepancy when a single parametrisation is applied, like the 



Phillips relation |Phillips! Il993h that relates a light-curve 
width to peak luminosi ty, due to in trinsically different explo- 
sions. Also, as noted in IWandl2005l . the circumstellar dust or 
a surrounding gaseous envelope may substantially affect the 
extinction properties compared to the interstellar dust. In 
other words, it is not wise to calibrate these SNe as standard 
candles; their reddenings will be more likely overestimated, 
which further implies that these objects are overcorrected to 
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a higher luminosity, that could bias the cosmological results. 
Moreover, a certain number of SNe fall within the rough 
range —0.1 < c < 0.1. Interestingly, imposing this cut-off 
criterion in Fig. 1(a) the correlation between scale 15 co- 



efficient and colour appears to be weak and the existence 
of distinct SN groups is no longer maintained. Nevertheless, 
one could select a subset of "more standardisable" SNe la 
based on the value of the largest wavelet scale coefficient. 

3.1 Intrinsic SN colour and host galaxy reddening 

As previously seen in Fig. [TJ the scale 15 coefficient is corre- 
lated with both the SALT2 colour and long-baseline B — I 
colour in a quite similar manner. Moreover, a strong rela- 
tionship between these latter two is shown in Fig. 2(a) The 
best fit line corresponds to: 

c ( salt2) = (0.114 ± 0.007) + (0.445 ± 0.012) x (B — I). (3) 



In Fig. 1(b) several SNe, namely 1990O, 1992A, 1994D, 



1996X, 1998aq, 1999aw and 2000dk, that are thought to suf- 
fer negligible extinction in their host galaxies (as compiled 
from the literature, Ay host ?S 0.04) are shown as crossed 
open circles. Putting all of these together, 

(B — I)max — E{B — I)mW = (B — /)max +E(B — /)host, 



(4) 



one expects that their colours corrected for Galactic extinc- 
tion, (B — I)™"', are equal to the unreddened intrinsic 
colours at maximum, (B — /)J£ax- The best fit that provides 
the intrinsic (B — I) colour relation for normal SNe is given 
by: 

(B - /)|'f ax = (0.013 ± 0.079) + (0.148 ± 0.029) x wltcoeffi 5 . 

(5) 

The extinction affects colour estimation and is a source 
of systematic uncertainty. The SALT2 light-curve fitter is 
not meant to be able to disentangle the potential extinction 
by dust in the host galaxies from intrinsic colour variation. 
In the case of nearby SNe the SALT2 colour law should be 
consistent with observed law f or extinction by dust in the 
Milky Way (extinction law by ICardelli. Clayton fc Mathisl 
1 19891 ). but SN data point to a diffe rent colour- luminosity re- 
lation (see e.g. lAstier et al. 1120061 '). It is interesting to note 
that in the case of SNe 1999aa and 1997cn, two SNe that be- 
long to other SN subsample discussed above and that do not 
suffer strong host galaxy extinction (Av = 0.025±0.020, i.e. 
0.025 ± 0.021 respec tively, according to the MLCS2k2 fits in 
Irlicken et al. |[2009l ). the relation given by Eqn. [5] is clearly 
unsati sfactory. Neverthele ss, this is in agreement with find- 
ings in I Wang et al~ll2009l . who reported two SN subsamples 
that favour different reddening laws. A possible explanation 
for this might be the existence of two SN populations with 
different intrinsic colours. 

Furthermore, expression in Eqn. [4] and the obtained fits 
allow us to estimate the host galaxy reddening, implicitly 
assuming that E(B-I) = 2.35E(B-V). Adopting Eqn.[5]as 
the intrinsic colour relation for (all) SNe, the value of E(B — 
V r )hoat can be directly obtained substituting expressions in 
Eqns. [3] and [5] in Eqn. H 

E{B - V) host = (0.069 ± 0.009) + (0.584 ± 0.028) x c (SALT2) , 

(6) 

including all SNe in the fit (shown as dashed line in 
Fig. 2(b) I. On the other hand, excluding SNe that do not 
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Figure 3. Comparison between wavelet coefficients at scale 15 
calculated on the high-z SNLS3 spectra for different wavelength 
intervals. The reference interval is 3400 — 7000 A. 



exhibit same properties as "normal" SNe (represented by 
filled symbols in Fig. 2(b) I from the fit and keeping only 
them in the second case, one gets: 



E{B - V)host = (0.035 ± 0.006) + (0.569 ± 0.018) x c {SALT2) , 



E(B-V) hoBt = (0.151 ±0.007) + (0.553 ± 0.023) xc {SALT2) , 
respectively. 

Again, there is a significant difference (~0.12 mag) in 
the reddening estimates between these two subsamples, im- 
plying that the two SN groups favour different reddening 
laws. However, this discrepancy may be due to assumption 
of the single intrinsic colour relation. Accounting for differ- 
ent intrinsic colour relations for the two SN subsamples may 
be a possible explanation for this. The main weakness of this 
analysis is the very limited number of normal SNe that en- 
ter the fit in Eqn. [5] also the lack of SN data that suffer 
no (or negligible) host extinction from other subsample that 
could provide an analogous intrinsic relation for intrinsically 
redder SNe. 



3.2 High-z sample 

An important question that needs to be asked is whether 
high-z SNe follow the same correlations as those in the pre- 
vious subsection. As an illustration, from the available set of 
high-z spectra from Super Nova Legacy Survey 3-year sam- 
ple (SNLS3), published in lBalland et al. 1120091 the wavelet 
analysis is applied to 94 SNe. Fig. [3] shows that the coef- 
ficient from the largest scale in the DTWT, measured on 
the high-z sample, is relatively stable to changes in upper 
wavelength limit, probably because the majority of high-z 
spectra do not reach 7000 A. Probing different wavelength 
intervals and identifying outliers, we find a high degree of 
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Figure 4. Comparison between wavelet coefficients associated 
with the largest scale calculated on the SALT2 spectral templates 
for different epochs and wavelength intervals after application of 
some amounts of reddening. 



agreement when comparing them. However, the wavelet co- 
efficient at largest scale increase^ if the lower range limit 
shifts towards larger wavelengths compared to the reference 
one of 3400 A, which might be an important issue when 
studying nearby SNe. 

3.3 Stability of the largest wavelet coefficient 

The results in Fig.[3]are compatible with tests based on mea- 
surements on the SALT2 spectral templates, as illustrated 
in Fig. 3] The wavelet coefficient from the coarsest scale is 
tested both for stability on chosen wavelength interval and 
epoch relative to B-band maximum. Although the disper- 
sion seen when comparing different measurements can be 
up to 50 per cent near maximum, the coefficient value is 
relatively invariant across epoch over the range [—10, +10] 
days. Furthermore, to include the dust effects, the spec- 
tral templates in the different wavelength intervals were 
reddened using the Milky Way (M W)-like extinction curve 
l|Cardelli. Clayton fc Mathisll 19891 ) for different values of Rv 
and E(B — V). It has been shown that the effect of reddening 
does not strongly affect the equivalen t width of a spectral 
feature, especially if it is narrow (e.g. iBronder et aLll2008l ; 
lArseniievic et a.lT1l2008h . Similarly, despite the overall spec- 
trum is somewhat warped, the relative error between dif- 
ferent cases becomes significant mainly for large extinctions 
(Fig. |4]). Furthermore, it can be seen that the dependence 
of the largest scale wavelet coefficient on extinction effects 
may be adjusted using specific wavelength intervals in the 
wavelet decomposition. 

To test the effects of reddening, a subset of the "normal" 
SNe la was randomly chosen and some amount of reddening 

2 and opposite for positive values 



was applied to them. Indeed, after obtaining new coefficients 
from the fits, the differences are sufficiently large that may 
indicate the existence of two populations of SNe. However, 
as we have seen in Fig. [4] using the right wavelength inter- 
val, these differences can be controlled, i.e. minimised. For 
instance, choosing the intervals such as [3700 — 6000 A] or 
[4000 - 6000 A], the maximum relative error reduces to a 
few per cent, even for large extinction values that have been 
added. 

The previously described method was performed on the 
same SN spectra, using the interval [4000 - 6000 A] (see 
Fig-©- This choice guarantees the stability of the calculated 
values of the wavelet coefficient to the reddening uncertain- 
ties or instrumental distortions. 

A robust regression analysis was performed to deter- 
mine the best fit to the data in Fig. |5(a)| 

c ( salt2) = (2.36 ± 0.12) + (0.73 ± 0.04) x wltcoeff 15 , (7) 

where the asterisk denotes the coefficient obtained in the 
interval [4000 - 6000 A]. 

Moreover, in Fig. |5(b) it is demonstrated how to per- 
form a test that identifies SNe la with anomalous extinction 
or strong peculiarity. For instance, apart from four outliers 
in Fig. [T(a)1 (SNe 1986G, 1999cl, 2003cg and 2006X), pe- 
culiar supernova 2008ha (extremely faint, 2002cx-like SN) 
strongly deviates from the m ain trend despite ne gligible host 
galaxy extinction reported in lFolev et al. |[2009l . A few high- 
z SNLS SNe whose spectra cover wavelengths up to 6700 A 
are added to the plot (given as filled diamonds) showing an 
excellent agreement with nearby SNe. 

Finally, different amounts (0 < Av < 4, with the step 
AAy = 0.4) of the MW-like extinction are applied to the 
SALT2 spectral template at maximum and obtained coeffi- 



cients are plotted in Fig. 5(b) The data points follow the 



usual SN trend, even for the extinctions as large as 2.5 mags. 
Nevertheless, the figure clearly shows the two regimes, a lin- 
ear (Av < 2.5), while there is an increase to the non-linear 
regime for heavily obscured SNe with Av > 2.5. In addi- 
tion, strongly extinguished SNe such as 1999cl, 2003cg and 
2006X, are consistent with this higher power law. After test- 
ing different extinction laws, we conclude that there is de- 
generacy between different models in the linear region; the 
differences mainly occur at higher values of extinction, i.e in 
the non-linear regime. The dependence of this combination 



of wavelet coefficients shown in Fig. 5(b) can be certainly 



used to estimate the amount of extinction towards an indi- 
vidual supernova. 



4 DISCUSSION 

We are driven by need to grasp whether the existence of 
spectroscopi c dichotomy, or e ven a quadrotomy, recently 
discu ssed in IWang et al. 1 l2009l and iBranch. Dang fc Baronl 
I2009I respectively, implies either different ext inction laws 
or different c olour evolution. As suggested in IWang||2005l ; 
iGoobari l200sl . the interaction between the SN ejecta and 
the surrounding circumstellar shell with normal dust grains 
could explain the reported values of Rv that are lower than 
the standard Galactic value of 3.1. Furthermore, the SN 
spectra at early epochs are sensitive to the material in which 
the supernovae are embedded. Emitted SN light is therefore 
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Figure 5. (a) Wavelet coefficient from the largest scale measured on the wavelet interval [4000 — 6000 A] versus SALT2 colour, (b) 
Comparison between the coefficients for different wavelength intervals. A linear fit (dashed line ) can be used to test SNe la for anomalous 
extinction and/or peculiarity. 



reflected by the surrounding dust clouds and reaches the ob- 
server at different times from different directions, as a result 
of the interaction that, generally speaking, depends on the 
location of the dust clouds, the density distribution of dust 
grains and/or the scattering properties of dust. Moreover, if 
the dust optical depth is low, one may assume a simple scat- 
tering scenario in which a photon escapes after scattering 
off dust particle. It becomes more difficult when thick dust 
shells and complex dust mixtures are con sidered, where mul- 
tiple scatterings usu ally occur. Following IPatat et al~ll2007l : 
iBlondin et"aLl2009l . the presence of high- velocity features in 
spectra and time- variable Na ID absorption feature is good 
signature of influence of effects of circumstellar dust on the 
SN luminosity (typical examples are SNe 1999cl and 2006X). 
Besides, SN 2006X suffers significant r eddening from the 
interstellar medium (|Wang et al. I I200H ). However, not all 
highly reddened SNe la display variable Na ID features: SN 
2003cg is heavily extinguished, very large Na ID equivalent 
width values (compared to other SNe) were measured for 
this SN, but the feature remains constant over time. The 
previously described scenario can cause an effect of the light 
echo around supernova, already seen, for inst ance, for SNe 
1991T 1995E, 19 98bu, 2000cx, 2002ic, 2006X (|Patat et al. I 
120071 ; I Wane! 120051 ). Similarly, the dust shells that surround 
(not necessarily highly extinguished) SNe la, are reported 
also for SN 2005g£], then 2007gi and 20071e. Indeed, all 
of these mentioned SNe are leaned towards the second fit 
in Fig. 1(a) This latter s upernova, recently discovered SN 
Simon et al. 120091 ). however, exhibits variable Na ID 



20071c 



consider the interaction between the ejecta and the circum- 
stellar material as an alternative manner for a supernova to 
generate extra luminosity, rather than any intrinsic property 
of the SN. Moreover , a supernova with extremely low lumi- 
nosity, SN 2OO8h£0 (|Folev et al. 1120091). that is spectr oscop- 
ically similar to peculiar SN 2002cx ( Li et alTI 120031 ). is in 



line but, curiously, is low extinguished, which further re- 
flects its intrinsic property. In addition, lYuan et al. 2010 



an excellent agreement with the red solid line in Fig. 1(a)! 
altho ugh it is an outlier in Fig. \E\ As noted in iLi et al.H 
I2003L SN 2002cx is a link between the extremes of peculiar 
SNe la. Indeed, this SN exhibits the 1991T-like premaxi- 
mum spectra, but is subluminous, such as 1991bg-like SNe. 
Nevertheless, the wavelet coefficient from the largest scale 
reflects spectroscopic diversity of SNe la, pointing to a com- 
mon property of this trend, regardless if an SN has already 
been classified as 1991T-, 1991bg-, 1999aa- or 2000cx-like 
event . 

On the other hand, SNe such as 1999ee and 1989B 
are also found to show affinity for this group. These SNe 
exhibit a plateau in the photospheric velocity near B- 
maximum, while normally, a standard SN la has very 
steep velocity change from early phases towards maxi- 
mum. This plateau tendency (more precisely, a period 
of slowly declining velocities) coincides wit h the photo- 
sphere receding through the overdense shell jScalzo et al. I 
l20ld ; IQuimbv, Hoflich fc Wheeler] l2007h . The similar has 
been noticed for SNe 1991T, 1999aa, 2000cx (see e.g. 
iBenetti et al~ll2005l ; IQuimbv. Hoflich fc Wheeled \200H ), all 
of them having redder (B — V) colour at maximum. 
IQuimbv. Hoflich fc W heeler 2007 report nearly the same ve- 
locity evolution of the Si II line for SN 2005hj, but the avail- 
able spectra for this SN do not cover wavelengths A < 3885, 



3 The SALT2 fit for this unusual SN is relatively poor, but we 
will use the obtained estimates given in Table [JJ 



4 Given as an open circle with high colour value, c = 0.411 in 
Fig- Eg] 
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therefore the value of wavelet coefficient of interest is most 
likely overestimated (see Figs. [3] and [4]). 

Another intriguing result is that the high velocity gra- 
dient (HVG. lBenetti et al. lliooH l SNe, namely SNe f998de, 
f999gd, 2002bo, f984A, 200 2bf, 1 991M and 2004dt, given 
as light filled circles in Fig. |f (a)| and ordered by decreas- 
ing colour, are found to be transitional objects, linking the 
two main SN trends. A possible interpretation of this ob- 
served continuity between the groups could be accounted as 
an effect of mixed extinction that these SNe suffer, partly 
caused by circumstellar, but could also be due to inter- 
stellar dust. However, another HVG SN 1981B is found in 
the CSM group. In addition, a few more SNe are identi- 
fied very close to these HVG SNe in direction towards nor- 
mal SNe, namely: f998bp, f998bu, f999aw, 2000cf, 2001V, 
2001el. Again, these latter are characterised by whether the 
higher polarisation (e.g. SN 2001el), or close resemblance to 
1999aa-like or under(over)lu minous SNe. On the other hand, 
it has been found recently (|Maeda et al. 1 120101 ") that HVG 
and LVG SNe, despite the divergence in photospheric veloc- 
ity, do not have intrinsic differences. The diversity between 
these two supernova groups appears as a consequence of dif- 
ferent viewing angle from which an SN is observed. In other 
words, if one obtains a big enough SN sample, these differ- 
ences will average out and this issue will be of no concern 
for the cosmological use of SNe la. However, the redder SN 
population that satisfies relation in Eqn. [2l must be treated 
with caution because our findings point to different intrinsic 
colour relations for the two SN subsamples. 

Given the previous discussion, our analysis suggests the 
existence of a distinct class of SNe that is more likely gov- 
erned by different physical mechanism, i.e. considered to 
have a different progenitor system and/or explosion mech- 
anism than standard SNe la. Intere stingly, SNe that sup- 
port both single- (e.g. SN 2005gj, lAlderfng etaTI l200rj ; 
iPrieto et aL~ll2007l ) and double-degenerate (e.g. SN 2003fg 
and/or SN 2006gz, with detected envelope of unburned car- 
bon) scenario are found in the same group. Nevertheless, if 
one suspects that the circumstellar dust is present around an 
SN la, it is not advisable to assume a standard interstellar 
extinction law, but rather should be studied it by imple- 
menting a radiative transfer code. The effective extinction 
is supposed to be generally smaller when the scattering ef- 
fects are included into consideration, since the contribution 
of scattered photons makes that the total flux increases. The 
light that has been scattered off the circumstellar dust cloud 
aims to lower the effective total to selective extinction ratio 
in the optical range, compl etely opposi te from the effect it 
causes m the UV (see e.g. I Wand 120051 ). A possibility that 
SN colours correlate with intrinsic brightness of SNe la in 
such a way that this correlation resembles reddening by dust 
is not rejected as well, as has been indirectly shown in rela- 



tions from Fig. 1(a) and Eqn. [5] It is also implicitly assumed 



throughout that the SALT2 colour estimation is essentially 
perfect, despite the fact that its training procedure is based 
on a sample of SNe la. 

Summing up, around 30 per cent of nearby SNe in our 
sample are found to depart from the main SN trend in 



rections for normal and high velocity (HV) samples, which 
notably improved the accuracy of distance measurements. 
Furthermore, the same authors showed that the HV SNe 
prefer a lower value of the extinction ratio Rv ■ On the other 
hand, our measurements on the spectral template show that 
the largest wavelet scale coefficient is fairly stable under dif- 
ferent extinction laws, which suggests that this observed dis- 
crepancy has an intrinsic character. Indeed, the two SN sub- 
samples may have the same host galaxy reddening law if the 
intrinsic colour relations were properly accounted for. 

In addition, even the highly-reddened SNe are con- 
sistent with the MW-like dust extinction behaviour in 
Fig. |5(b)| Tests with different Rv values than the canonical 
Galactic value of 3.1 also suggest that anomalous behaviour 
is most likely due to high r eddening. Th e extinction law for 
the circumstellar dust from lGoobarll200"8l is also tested and it 
appears very much like the curve with the MW extinction in 
the linear part, while it is steeper in the non-linear regime, 
although for significantly smaller values of extinction (for 
example, for a fixed value of Rv = 1.8, a CS curve that fits 
SN 2003cg measurements has Av ~ 2.2, compared to the 
MW curve with Av ~ 2.5). In order to include a more sig- 
nificant high-2 SN sample into the analysis (compared with 
those few SNe that are found in Fig. 5(b) |, another pair of 



Fig. 1(a) Assuming that the similar fraction of high-z pop- 
ulation follows this tendency, a contamination of cosmologi- 
cal sample by thes e objects could be significant. Recall that 
I Wang et al.Tl2009l recently proposed two different Rv cor- 



the wavelength intervals should be studied, ones that corre- 
spond better to the rest-frame coverage of high-z SN spectra, 
for example between 3000 < A < 5000 A. 

Wavelet scales that are directly responsible for spec- 
tral features are free of any intrinsic color informa- 
tion in the input spectra, which justifies a recent at- 
tempt to quantify SN spectral features that was done by 
IWagers. Wang fc Asztaloal2010l . They found correlations be- 
tween various combinations of spectral features with stretch. 
This is surely an important question that should be further 
explored. 



5 CONCLUSIONS 

The presented technique is used to test SN diversity, i.e. 
how significantly the explosion of a supernova differs from 
widely accepted SN la progenitor scenarios, which may fur- 
ther lead to different intrinsic colour. For this purpose, the 
wavelet transform is used to decompose SN spectra into dif- 
ferent scales. The largest scale coefficient is found to cor- 
relate with the SALT2 colour parameter and long-baseline 
B—I colour corrected for Galactic reddening. Apart from the 
main trend with normal SNe la, another grouping is distin- 
guished; its members have intrinsically redder colours when 
compared to the normal ones and are recognised for vari- 
ous characteristics such as the observed light echo around 
SN, or in addition, a signature of interaction with circum- 
stellar material, and/or the existence of shell-structure in or 
surrounding the SN ejecta. Moreover, a few candidates of 
super-Chandrasekhar mass SNe show affinity for this group. 
However, the questions on details on circumstellar dust shell 
(or inner overdense shell), such as its mass or its geometry, 
are not tackled in this work. 

This study has shown that SALT2 colour parameter can 
be disentangled into intrinsic SN colour and the reddening 
due to dust extinction in the host galaxy. However, with 
the lack of SNe from intrinsically redder subsample that 
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suffer negligible host galaxy extinction, the expressions for 
intrinsic colour and host galaxy reddening have been found 
only for normal SNe. In order to prevent the overcorrected 
luminosities, mentioned redder SNe should be identified and 
calibrated in a different manner before being included in the 
cosmology fits. 

The wavelet coefficients that were measured on certain 
wavelength intervals are invariant to additional amounts of 
extinction that were applied to the spectra. A combination 
of two largest scale coefficients from different intervals can 
be further used to distinguish objects that exhibit strong 
peculiarities and/or to estimate the extinction value. 

The method for additional SN sorting presented in this 
paper is applicable regardless either of the type of wavelets 
used in the analysis or of the number of scales in the wavelet 
decomposition. Nevertheless, there is certainly room for im- 
provement, especially in estimating the largest wavelet scale 
coefficient when the available spectra lack the full wave- 
length coverage. Similarly, an additional correction could be 
applied to account for spectral epochs different from the 
reference one at maximum light. The implementation of the 
presented analysis by incorporating larger and homogeneous 
spectral and photometric data will permit a deeper insight 
into the extinction corrections, also the variety of the SN 
progenitor environments and explosion models and will be 
the topic of future work. 
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Table 1: List of SNe la used in our study. 



SN name 


Redshift 






Xl 




c 


X 7dof 




wltcoeffis 




wltcoeffis 


Refs 


1981B 


0.0072 


-0 


.472 ± 0. 


,138 





,048 ± 0.038 


1.032 


-4 


,428 ± 0.800 


-2 


.948 ± 0.624 


1 


1984A 


-0.0009 





,051 ±0 


,050 





,218 ±0.032 


26.152 


-1 


,622 ± 1.304 


-2 


.555 ±0.147 


4 


1986G 





.0018 


-2 


.226 ±0. 


.142 





,845 ± 0.042 


1.378 


3 


,128 ±0.639 


-1 


,536 ± 0.256 


1 


1989B 





.0024 


-0 


,913 ±0. 


.105 


0. 


,343 ± 0.028 


0.713 


-1 


,767 ± 1.565 


-3 


,181 ±0.589 


1 


1990N 





.0034 





,996 ± 0, 


.094 


0. 


,008 ± 0.025 


2.025 


-1. 


,680 ± 0.349 


-3 


,242 ± 0.692 


1 


1990O 





.0307 





,450 ± 0, 


,191 


-0. 


,053 ± 0.032 


0.296 


-2. 


,740 ± 0.500 


-3 


,037 ± 0.058 


10 


1991M 





.0072 


-1 


,082 ± 0, 


239 


0. 


,018 ± 0.089 


0.569 


-3 


,254 ± 0.800 


-3 


,103 ±0.596 


1 


1991T 





.0058 


1 


,410 ±0, 


,087 


0. 


,112 ± 0.025 


2.100 


-3 


,718 ±0.800 


-3 


,235 ±0.760 


1 


1991bg 





.0030 


-1 


,377 ±0, 


,057 


0. 


,290 ± 0.034 


31.701 





,118 ±0.800 


-2 


,466 ±0.614 


1 


1992A 





.0061 


-1 


,632 ±0, 


,051 


-0. 


,013 ± 0.024 


1.368 


-2 


,115 ± 0.326 


-3 


,136 ±0.514 


1 


1994D 





.0015 


-1 


,741 ± 0, 


,045 


-0 


,095 ± 0.024 


1.258 


-2. 


,893 ± 0.346 


-3 


,244 ±0.381 


1 


1994U 





.0048 





,000 ± 7, 


,465 


0. 


,297 ± 0.323 


3.341 


-2. 


,507 ± 0.800 


-3 


,447 ± 0.660 


14, 15 


1995E 





.0118 


-0 


,488 ± 0, 


122 





,723 ± 0.030 


0.542 


-0. 


,072 ± 0.443 


-3 


,223 ± 0.576 


15, this work 


1996X 





.0068 


-1 


Oil ±0, 


,095 


-0. 


,039 ± 0.028 


0.200 


-2 


,582 ± 0.277 


-3 


,179 ±0.529 


1 


1997br 





.0070 





,369 ± 0, 


,098 


0. 


,200 ± 0.027 


7.670 


-0. 


,873 ± 0.206 


-3 


736 ± 0.661 


1 


1997cn 





.0162 


-0 


,550 ±0, 


.116 


0. 


,777 ± 0.030 


23.903 


1. 


,007 ± 0.800 


-2 


,274 ± 0.467 


1 


1997do 





.0101 





,220 ±0, 


.155 


0. 


,070 ± 0.026 


2.567 


-2 


,020 ±0.551 


-3 


,024 ± 0.497 


2 


1997dt 





.0073 


-0 


,141 ±0, 


.290 


0. 


,574 ± 0.080 


0.346 





111 ±0.066 


-3 


,602 ± 0.558 


2 


1998V 





.0176 


-0 


,333 ± 0, 


,087 


-0. 


,004 ± 0.026 


1.734 


-1 


,951 ± 1.296 


-3 


,149 ± 0.673 


2 


1998aq 





.0037 


-0 


,353 ±0, 


036 


-0. 


,150 ±0.025 


1.398 


-2 


,848 ± 0.353 


-3 


,088 ± 0.672 


1, 2 


1998bp 





.0104 


-2 


,440 ± 0, 


,169 





,184 ±0.031 


1.108 


-1 


,618 ±0.435 


-3 


,246 ± 0.644 


2 


1998bu 





.0030 


-0 


,344 ± 0, 


,029 


0. 


,185 ±0.027 


3.116 


-1. 


,671 ± 0.458 


-3 


,266 ± 0.674 


1, 2 


1998de 





.0167 


-3 


,061 ±0, 


133 





,443 ± 0.029 


2.988 


-0 


,091 ± 0.684 


-2 


,341 ± 0.465 


2 


1998dh 





.0089 


-0 


,692 ± 0, 


,106 


0. 


,055 ± 0.027 


1.346 


-2. 


,062 ± 0.123 


-3 


,169 ±0.346 


2 


1998dm 





.0066 





,596 ± 0, 


.251 


0. 


,245 ± 0.031 


0.965 


-0 


,621 ± 1.160 


-3 


,585 ±0.427 


2 


1998ec 





.0199 


-0 


,188 ±0, 


323 


0. 


,134 ±0.057 


0.435 


-1 


,153 ±0.582 


-3 


,121 ±0.652 


2 


1998eg 





.0248 


-0 


,700 ± 0, 


316 


-0. 


,009 ± 0.028 


0.370 


-2 


,151 ±0.329 


-3 


,068 ±0.631 


2 


1998es 





.0106 


1 


Oil ±0, 


.119 





,036 ± 0.025 


0.553 


-4. 


,600 ± 1.523 


-3 


,144 ±0.674 


2 


1999X 





.0252 


-0 


,529 ±0, 


.498 


0. 


,010 ±0.050 


0.271 


-2. 


,293 ± 0.191 


-3 


,096 ± 0.524 


2 


1999aa 





.0144 


1 


,170 ±0, 


,093 


-0 


,079 ± 0.025 


0.803 


-4, 


,817 ± 1.470 


-3 


,486 ± 0.566 


1, 2 


1999ac 





.0095 





,206 ± 0, 


,064 


0. 


,029 ± 0.024 


3.220 


-2. 


,129 ± 1.388 


-3 


329 ±0.691 


1, 2 


1999aw 





.0380 


2 


,339 ± 0, 


166 


-0. 


,075 ± 0.029 


1.649 


-3. 


,425 ± 0.132 


-3 


,484 ± 0.033 


1,2 


1999by 





.0021 


-1 


,444 ± 0, 


063 


0. 


,360 ± 0.025 


17.604 


-2. 


,085 ±1.116 


-2 


,261 ±0.497 


1, 2 


1999cc 





0313 


-1 


,582 ±0, 


,199 


-0. 


,025 ± 0.031 


0.220 


-2. 


,826 ± 0.224 


-2 


,999 ± 0.627 


2 


1999cl 





.0076 


-0 


,753 ± 0, 


,167 





,984 ± 0.027 


1.581 


3. 


,622 ± 0.582 


-1 


,726 ± 0.303 


2 


1999dq 





0143 





,826 ±0, 


,062 


0. 


,033 ± 0.025 


0.867 


-2 


,193 ± 2.014 


-3 


,457 ± 0.682 


2 


1999ee 





.0113 





,722 ±0, 


039 





,231 ±0.023 


0.690 


-2 


,888 ±0.191 


-2 


,785 ± 0.576 


1 


1999ej 





.0137 


-1 


,400 ± 0, 


.217 


-0. 


,037 ± 0.038 


0.490 


-2 


,290 ± 0.272 


-3 


,206 ± 0.493 


2 


1999gd 





.0185 


-1 


,082 ± 0, 


126 


0. 


,372 ± 0.029 


1.829 


-0. 


,555 ± 0.070 


-3 


,134 ±0.629 


2 


1999gh 





.0077 


-2 


,329 ±0, 


.152 





,120 ± 0.034 


0.381 


-0. 


,828 ± 0.410 


-2 


,554 ± 0.478 


2 


1999gp 





.0267 


1 


,749 ± 0, 


.083 


-0. 


,004 ± 0.018 


1.082 


-2. 


,539 ± 0.936 


-3 


,068 ±0.619 


2 


2000E 





.0047 





,476 ± 0, 


,062 





,102 ± 0.024 


3.382 


-1. 


,839 ± 0.348 


-3 


,528 ± 0.668 


1 


2000cf 





.0364 


-0 


,959 ±0, 


,210 


-0. 


,034 ± 0.029 


0.467 




,245 ±0.101 


-2 


,967 ± 0.392 


2 


2000cn 


() 


.0235 


-2 


,247 ± 0, 


,209 


0. 


,077 ± 0.028 


0.652 


-1. 


,562 ± 0.630 


-3 


,311 ±0.114 


2 


2000cx 





.0079 


-1 


,001 ±0, 


,042 


-0. 


,088 ± 0.024 


5.622 


-4. 


,849 ± 0.207 


-3 


,430 ± 0.550 


1, 2 


2000dk 





.0174 


-1 


,871 ±0, 


,144 


-0. 


,032 ± 0.027 


0.454 


-2. 


,194 ± 0.265 


-3 


,335 ± 0.439 


2 


2000fa 





.0213 





,554 ±0, 


,187 


0. 


,035 ± 0.029 


0.249 


-2 


,292 ± 0.771 


-2 


,927 ±0.637 


2 


2001V 





0160 





,988 ± 0, 


,094 


-0. 


,010 ±0.020 


0.950 


-3. 


,091 ± 0.675 


-3 


,212 ±0.607 


2 


2001el 





.0039 


-0 


,245 ± 0, 


,039 


0. 


,097 ± 0.024 


3.477 


-2. 


,487 ± 0.826 


-2 


,830 ± 0.050 


1 


2002bf 





.0242 


-0 


,298 ± 0, 


.145 


0. 


,144 ±0.034 


1.087 


-2. 


,268 ± 0.800 


-3 


,170 ±0.493 


1 


2002bo 





.0042 


-0 


,471 ± 0, 


.038 


0. 


,339 ±0.018 


1.425 


-0. 


,909 ± 0.951 


-3 


,157 ±0.615 


1 


2002cx 





.0240 





,351 ±0, 


136 





,308 ± 0.031 


6.960 


-2. 


,620 ± 0.400 


-3 


,082 ± 0.565 


1 


2002dj 





.0094 





,028 ±0, 


,146 





,060 ± 0.025 


0.934 


-3 


,781 ± 0.568 


-2 


,831 ±0.631 


13 


2002er 





.0086 


-0 


,764 ± 0, 


,057 


0. 


,100 ± 0.024 


0.929 


-1 


,112 ± 1.054 


-3 


,423 ± 0.620 


1 


2003cg 





.0041 


-0 


,802 ±0, 


039 


1. 


,053 ±0.019 


4.988 


5. 


Oil ±0.312 


-1 


,249 ± 0.204 


1 


2003du 





.0064 





,294 ± 0, 


039 


-0. 


,102 ± 0.017 


1.342 


-2. 


,686 ± 0.187 


-3 


,058 ± 0.437 


1 


2003fg 





.2437 


1 


,134 ±0, 


,306 


0. 


,043 ± 0.030 


10.955 


-4. 


,459 ± 0.500 


-3 


,824 ± 0.638 


12 


2004S 





.0091 


-0 


,217 ±0, 


,087 


0. 


,022 ± 0.026 


0.933 


-2. 


,612 ± 0.800 


-2 


,912 ±0.481 


1 



continued on next page 
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Table [T] continued 



SN name 


Redshift 




Xl 




c 


x7dof 


wltcoeffis 




wltcoeffi 5 


Refs 


2004dt 


0.0197 


-0.885 ±0 


048 


-0 


027 ± 0.024 


3.012 


-3.567 ± 1.113 


-3 


090 ± 0.594 


1, 3 


2004eo 


0.0157 


-1.044 ± 


036 





039 ± 0.024 


1.098 


-1.982 ±0.136 


-3 


549 ±0.155 


1 


2005am 


0.0079 


-2.034 ± 


094 





133 ± 0.040 


1.642 


-1.096 ±0.800 


-.3 


755 ± 0.593 


1 


2005bl 


0.0251 


-2.185 ±0 


100 





513 ±0.027 


7.728 


-0.667 ± 0.653 


-1 


465 ± 0.342 


11 


2005cf 


0.0065 


-0.134 ±0 


054 


-0 


034 ± 0.024 


1.119 


-2.356 ±0.149 


-3 


235 ± 0.463 


1 


2005gj 


0.0593 


5.000 ± 


015 





178 ± 0.034 


9.610 


-3.464 ± 0.800 


-2 


941 ± 0.587 


19,20 


2005hj 


0.0580 


2.279 ± 


408 


-0 


077 ± 0.032 


0.778 


-1.703 ±0.967 


-3 


679 ± 0.094 


1 


2005hk 


0.0130 


0.044 ± 


058 





250 ± 0.027 


32.085 


-2.808 ± 0.390 


-3 


537 ± 0.670 


1 


2006D 


0.0097 


1.341 ±0 


044 





019 ±0.023 


1.288 


-3.645 ± 0.205 


-3 


200 ±0.164 


16 


2006X 


0.0063 


-0.258 ±0 


040 


1 


209 ± 0.025 


2.411 


4.434 ± 0.326 


-1 


667 ± 0.253 


5,6 


2006gz 


0.0280 


2.228 ±0 


120 


-0 


027 ± 0.025 


1.018 


-3.898 ± 0.347 
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